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We investigate the role of electron correlation effects in the electric dipole polarizabilities of the 
ground states of the alkaline earth and ytterbium atoms by employing the relativistic coupled- 
cluster (RCC) theory. These effects are incorporated via the residual Coulomb interaction to all 
orders in the RCC singles and doubles approximation. The perturbed wavefunctions used in the 
calculations of the polarizabilities are obtained by directly solving the first order perturbed RCC 
equations, thereby avoiding the sum-over-states approach. Our results are compared with other 
calculations and available experimental data. 
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I. INTRODUCTION 

A knowledge of electric dipole polarizabilities is neces- 
sary in many areas of physics and chemistry. In partic- 
ular it is required in studies of collisions involving atoms 
and molecules In recent years, the pre-eminent role 
of polarizabilities in the determination of inter atomic in- 
teractions has assumed special significance in the context 
of research on ultra-cold atoms [2] ■ 

Calculations of atomic polarizabilities have come a long 
way since the classic work of Dalgarno and Lewis [3| . Fol- 
lowing a series of calculations using the coupled Hartree- 
Fock method (see for example, |4j), a number of state- 
of-the art methods including the coupled-cluster (CC) 
theory have been used to calculate atomic and molecular 
polarizabilities [!, El 0] ■ A few calculations of the polar- 
izabilities of heavy atomic systems have been performed 
in the past few years using the linearized Q as well as 
the non-linearized @ relativistic CC theory. These are 
based on approaches that sum over a set of intermediate 
states. In contrast, we have obtained the first order per- 
turbed wavefunction by solving the first order perturbed 
CC equation and used it to obtain the dipole polariz- 
abilities of the closed-shell alkaline earth atoms and yt- 
terbium (Yb). Information on these quantities is useful 
for the frequency standards experiments that have been 
proposed for Mg [uj, Ca dH, Sr [l| and Yb [3 as well 
as the search for the electric dipole moment in Yb [3] • 

There has been considerable interest in accurate cal- 
culations of the dipole polarizabilities of alkaline earth 
atoms and Yb. Sadlej et al have calculated these quan- 
tities for Ca, Sr and Ba using a quasi relativistic ap- 
proach but their treatment of correlation is rigorous [la ] . 
Their calculations have been carried out at the finite or- 
der many-body perturbation theory and CC levels. Por- 



sev and Derevianko have performed calculations of the 
dipole polarizabilities of Mg, Ca, Sr and Yb by a hybrid 
approach involving many-body perturbation theory and 
the configuration interaction method 10]. A calculation 
of the dipole polarizability of Yb based on the time de- 
pendent density functional theory (TDDFT) has been 
reported recently [l6| . 



II. THEORY AND METHOD OF 
CALCULATIONS 



In a DC electric field E = £z, the energy shift AE 
of the ground state |<I/( ' (7, J , Mo) > with the parity 
eigenvalue 7 and angular momentum Jo and its azimuthal 
value Mo is given by 



AE = --a£ 2 , 
2 



(2.1) 



where a is the static polarizability and can be defined as 
9 \- \(^ Hl,Jo,M )\D z \¥°)(Y,J I ,M I ))\\ ^ 



where D z is the zth component of the electric dipole op- 
erator depending upon the applied field E, subscript / 
denotes the intermediate states and Jj, Mj are their an- 
gular momentum quantum numbers, 7 and 7' are parity 
quantum numbers for states of opposite parity and Eq 
and Ej are the energy values of the ground state and the 
corresponding intermediate states /, respectively. 

In a more explicit form, the above expression can be 
written as 
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(^(°)(7,Jo,Mo)|^|vI/( )(7 , ,Jj,M J ))(vl/( )(y,J J ,M J )|J z |vl/( )( 7 ,Jo,M )) 

So — Ej 



(vl/( )( 7; Jo ; M )|J z |^( )(Y,J J ,M J ))(vl/( )( 7 ',J J; Af J )|^|vl/( )( 7 ,Jo,M )) 
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where we define (7', J , M )) as the first order 

correction to the original unperturbed wavefunction 
>I>( )( 7 , Jo, Mq)) due to the operator D z and hence the 
total wavefunction in the presence of an external DC elec- 
tric field is given by 

|*(7,J 0j M )> = |* (0) (7,^o, M )) + |*«( 7 ',J , M )). 

(2.4) 



Our aim in this work is to obtain the exact wavefunc- 
tion *( 7 , J , M )) by calculating both * (0) (7, Jo, M )) 
and l^/*- 1 - 1 (7, Jo, Mo)} wavefunctions using an approach 
which can rigorously incorporate the relativistic and cor- 
relation effects . In other words, we would like to ob- 
tain |* (1) (7', Jo, M )) for the Dirac-Coulomb (DC) wave- 
function |'I'( )(7, Jo, Mq)) as the solution of the following 
equation 



(Ho° C) 



4°))|*W( 7 ', Jo, M )) = (£ -H int )|*(°)(7, Jo, Mq)), 



(2.5) 



r 



where H DC ' ) and Hj nt correspond to the DC Hamiltonian 



and interaction due to D z , respectively. 

Using the CC ansatz, we express 
|*( 7 ,J ,M )) by [jT 



explicitly, we can write the unperturbed and perturbed 
wavefunctions in terms of the RCC excitation operators 



wavefunctions as 



V± °" ^ ^ |* (0) (7, Jo, Mo)) = e T( > (7, Jo,M,)), (2.8) 

|*(7,J ,Mo)} = e T |$ (7,Jo,M )), (2.6) |*«( 7 , J 0) M )) = e r(0) T« |* ( 7> J , M )), (2.9) 



e T(0) |$ (7, Jo, M )), 



where $0(7, Jo, Mo)} are the Dirac-Fock (DF) wavefunc- 
tions determined using the mean-field approximation and 
T are the electron excitation operators from the corre- 
sponding DF states. 

To obtain both the unperturbed and perturbed wave- 
functions of Eq. (12. 4|) separately, we divide the excitation 
operators T as 



rp _ J>(0) _|_ ^y(l) 



(2.7) 



where and are the all order excitation operator 
of the relativistic coupled-cluster (RCC) method and its 
first order correction arising in the presence of D z , respec- 
tively. A represents the perturbation parameter . More 



where the exponential function of reduces to the 
linear term since we have considered only one order of 
D z operator in Eq. (2.5). 

In the present work, we have considered all possible 
single and double excitations (CCSD approach) in the 
calculations. We obtain first the unperturbed ampli- 
tudes by solving the usual (R)CC equations, then these 
amplitudes are used to determine the T^ 1 -* amplitudes. 
The corresponding equations to solve these amplitudes 
are given by 



(<frg|Hff c) |<I>o) = 
($S|H^ C) T«|$ ) 



~<$*|H int |$ ) 



(2.10) 
(2.11) 
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where the subscript N represents normal order form of 
the operators and we have defined H = e _T<0) He T< ' = 
(He T< )co n , with the subscript 'con' representing con- 
nected terms. 

The polarizabilities are determined by evaluating the 
following expression 

(4> |e Tt De T \<£ n ) 
<$ |e Tt e T |$ ) 

= (^oKrW'DW + pwrW)!^ > 

<$ |e T(0)t e T(0, |$o) ' 



where for computational simplicity we define D(°) = 

e T< °' De T< °'. We compute these terms after expressing 
them as effective one-body and two-body terms using the 
generalized Wick's theorem [181 ]. 



III. RESULTS AND DISCUSSIONS 

In Table HI we present our results of electric dipole 
polarizabilities and compare with those available in the 
literature. The error bars in the experimental results are 
large for all the systems and our results lie within them. 
The results of the calculations that are given in the above 
table are obtained using a variety of many-body theories. 
Our results are in reasonable agreement with them except 
for a few cases where they differ by more than 5%. 

Mitroy and Bromley [lj| have used oscillator strengths 
from a semiempirical approach to obtain these quanti- 
ties. Some results based on an ab initio method that 
combines the configuration interaction (CI) method and 
man y-b ody perturbation theory (MBPT) are available 
[if! l2dl . [HI- I n these calculations, the valence elec- 
trons correlation effects are evaluated by the CI method 
whereas the core electrons correlation effects are calcu- 
lated using the MBPT method. However, these calcula- 
tions consider the core-polarization effects in the frame- 
work of finite order MBPT, while our CC method takes 
them into account to all orders. Lim and Schwerdt- 
feger had employed the scalar-relativistic Douglas-Kroll 
Hamiltonian to determine some of these quantities [22| . 
In these calculations, they had shown the importance 
of the relativistic effects. There are also a few calcula- 
tions available for the Yb polarizabilities using the CC 
method [23|, [24| , where the atomic orbitals are evaluated 
using the molecular symmetries. Recently TDDFT was 
used for calculating the same quantities (l6j |. but this 
method treats exchange and correlation effects via local 
potentials. The main difference between the methods on 
which all these calculations are based and ours is that we 
calculate both the unperturbed and the first order per- 
turbed wavefunctions using a RCC approach that implic- 
itly takes into consideration all the intermediate states. 

To emphasize the importance of correlation effects in 
these calculations, we present the DF and the leading 



TABLE I: Static dipole polarizabilities in divalent atoms: Be, 
Mg, Ca, Sr, Ba and Yb (in au). 



Atoms 


Expts 


Others 


This work 


Be 




37.755", 37.69 b , 37.9 C 


37.80 


Mg 


71.5(3.1) d 


71.35\ 72.0 C , 71.3(7) e 


73.41 


Ca 


169(17) d 


159.4\ 152.7 C , 157.1(1.3) e 
152 / , 158.0 s , 


154.58 


Sr 


186(15) d 


201.2\ 193. 2 C , 197.2(2) e 
190 / , 198.9 s , 


199.71 


Ba 


268(22)'' 


264', 273 / , 273.9 s 


268.19 


Yb 


142(36) J 


111.3(5) e , 141.7 fe , 157.30' 


144.59 


Reference: f2| | 
b Reference: IS ] 
Reference: 2£ 
Reference: 27] 
e Reference: 1C 
■^Reference: 15 1 
Reference: 22 1 
^Reference: 2£ j 
'Reference: 21 
-'Reference: 2E 
k Reference: 2c 
'Reference: 16] 



TABLE II: Contributions from DF and important perturbed 
CC terms (in au) for the dipole polarizabilities. 



Atoms 


DF 


{Dt[ 1] + cc) -DF 


DT 2 (1) + cc 


Norm 


Be 


45.82 


-7.94 


-0.09 


0.02 


Mg 


82.44 


-8.77 


-0.21 


0.03 


Ca 


184.14 


-29.23 


-0.07 


-0.26 


Sr 


234.41 


-34.46 


-0.17 


-0.08 


Ba 


328.32 


-61.18 


0.09 


0.81 


Yb 


183.32 


-39.86 


0.032 


1.10 



RCC contributions in Table HT1 for the electric dipole po- 
larizabilities. For all the cases that have been considered, 
the DF results are larger than the total results. From 
the individual RCC contributions, we find that only the 
terms arising from DT^ and its conjugate (cc) are sig- 
nificant. Given that these terms include the DF, leading 
core-polarization and other important correlation effects 
to all orders, it is not surprising that they should col- 
lectively make up the largest contribution. In Fig [1] we 
give the breakdown of DT± in terms of the DF, ran- 
dom phase approximation (RPA) and other diagrams. 
All the above mentioned calculations where these terms 



4 



TABLE III: Contributions from important nsmp 3 Pi and 
nsmp 1 Pi states in the form of single particle orbitals through 
the (DT[^ + cc) CC terms to the dipole polarizabilities. 



Atoms ns mp DF DT^ + cc 



have been evaluated are based on finite order MBPT, 
but in our present method, we treat them to all orders 



in the residual Coulomb interaction. The 3 Pi and 1 Pi 
configurations are of crucial importance for the leading 
term of DT^ and its conjugate. In Table Hill we present 
the contributions arising from different combinations of 
single particle orbitals. 




(vii) (viii) fa) 



FIG. 1: Breakdown of the DT^ CC diagram in terms of 
lower order MBPT and RPA diagrams that contribute sig- 
nificantly to the polarizability calculations. Here, D and Vjv 
represent the dipole and normal order Coulomb interaction 
operators which are shown as single dotted and dashed lines, 
respectively. 



IV. CONCLUSION 

We have carried out calculations of electric dipole po- 
larizabilities for many alkaline earth atoms and ytter- 
bium using the RCC method and highlighted the impor- 
tance of the correlation effects. The novelty in these cal- 
culations is that it avoids the sum-over-states approach 
in determining polarizabilities and implicitly considers all 
the intermediate states by solving the perturbed RCC 
wavefunction to first order in the dipole and all orders in 
the residual Coulomb interaction. 
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